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Abstract 

We describe a simple formalism for generating classes of quantum 
circuits that are classically efficiently simulatable and show that the 
efficient simulation of Clifford circuits (Gottesman-Knill theorem) and 
of matchgate circuits (Valiant's theorem) appear as two special cases. 
Viewing these simulatable classes as subsets of the space of all quan- 
tum computations, we may consider minimal extensions that suffice to 
regain full quantum computational power, which provides an approach 
to exploring the efficacy of quantum over classical computation. 



1 Introduction 

The characterisation of the possibilities and limitations of quantum compu- 
tational power is one of the most interesting issues in quantum information 
science. All of the early and best known quantum algorithms [Ij that ex- 
hibit an exponential time speed-up over any known classical algorithm for 
the task, utilize properties of the quantum Fourier transform modulo A'^. 
One may then develop generalisations of these insights, studying Fourier 
transforms over further abelian and non-abelian groups and invent associ- 
ated computational tasks such as the hidden subgroup problem and various 
kinds of hidden shift problems. Around the years of 1997 and 1998 Thomas 
Beth, with memorable characteristic exuberance, was one of the earliest 
workers in the subject to recognise the potential possibilities of the abstract 
formalism of Fourier transforms for novel quantum algorithms, and take up 
this line of development which has now become an important cornerstone 
in our understanding. 

Despite this seminal development it is probably fair to say that apart 
from the Fourier transform formalism, no other similarly fruitful quantum 
algorithmic primitive for exponential speed-up has been identified. This mo- 
tivates a need for alternative approaches to exploring the efficacy of quantum 
vs. classical algorithms. One interesting such approach is the identification 
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and study of classes of quantum computations that are classically efficiently 
simulatable i.e. processes which although quantum, do not offer computa- 
tional benefit. Indeed the relation of classical to quantum computation that 
emerges is intriguingly rich and multi-faceted - (sub-) classical computation 
can be embedded into quantum computation in many inequivalent ways. 
Given any such class of simulatable quantum computations we may ask: 
what kind of added (minimal) ingredient suffices to restore full quantum 
computational power? In a sense, any such ingredient may be viewed as an 
"essence" of quantum computational power, albeit relative to a given sub- 
strate of simulatable processes. In this talk we will outline a formalism for 
providing simulatable classes of quantum circuits and discuss two examples 
- the Gottesman-Knill theorem for Clifford circuits and Valiant's theorem 
for simulation of matchgate circuits. These examples will show that the 
added ingredient above can be strikingly trivial, especially if thought of as a 
competitor to the oft-quoted blanket attribution of quantum computational 
power to the enigmatic phenomenon of entanglement. 

2 Classically simulatable quantum computations 

We focus on comparing and contrasting two theorems of classical simulation 
viz. the Gottesman-Knill theorem for Clifford circuits [U [9] and Valiant's 
theorem ^ [2] for simulation of matchgate circuits. At first sight these 
appear to be very different in their content and provenance but we will 
outline a proof method that reveals a formal similarity between the two 
results. 

The Gottesman-Knill (GK) theorem arose out of the development of 
the so-called stabiliser formalism for the theory of quantum error correction 
[1]. Let H denote the 1-qubit Hadamard gate, P the 1-qubit phase gate 
P = diag(l, i) and CZ the 2-qubit controlled— Z gate CZ = diag(l, 1, 1, —1). 
These gates and arbitrary circuits of them on n qubits are called Clifford 
operations on n qubits. Our adopted version (slightly modified from the 
original, c.f. also [3j) of the GK theorem is the following. 

Theorem 1 Consider any uniform (hence poly sized) quantum circuit fam- 
ily comprising the gates H,P and CZ (i.e. a Clifford circuit) such that: 
(i) the input state is any product state; 

(a) the output is a final Z measurement on any single qubit line. 
Then the output may be classically efficiently simulated. 

More formally our notion of efficient classical simulation is the following: 
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given a description of the circuit on n qubit lines, the output probabihties 
may be classicahy computed to k digits in poly(n, k) time. 

Next we introduce the notion of "matchgate" and Vahant's classical sim- 
ulation theorem [4], which arose originally from considerations of counting 
perfect matchings in graphs. 

A matchgate [U [2] is defined to be any 2-qubit gate G{A, B) of the form 
(in the computational basis): 



where A and B are both in SU{2) or both in U{'2) with the same determi- 
nant. Thus the action of G{A, B) amounts to A acting in the even parity 
subspace (spanned by |00) and |11)) and B acting in the odd parity subspace 
(spanned by |01) and |10)). 

Our version of Valiant's theorem (again slightly different from the origi- 
nal version) is the following. 

Theorem 2 Consider any uniform (hence poly-sized) quantum circuit fam- 
ily comprising only G(yl, B) gates such that: 
(i) the G{A,B) gates act on nearest neighbour (n.n.) lines only; 
(a) the input state is any product state; 

(Hi) the output is a final measurement in the computational basis on any 
single line. 

Then the output may be classically efficiently simulated. 

Let us now return to the GK theorem and its proof ingredients. The es- 
sential property of the class of gates used, i.e. Clifford gates, is the following 
[9]: if C is any n-qubit Clifford operation and Pi . . . ® P„ is any product 
of Pauli matrices (i.e. Pi = I,X,Y or Z for each i) then the conjugate 
C^(Pi . . . (g) Pn)C = P( . . . is again a product of Pauli operations. 
Stated more formally, if Vn is the group generated by all such Pauli products 
on n qubits then the n-qubit Clifford group is the normaliser of Vn in the 
unitary group ?7(2"). 

A standard proof (c.f. [Jj) of the GK theorem (with a computational 
basis input) proceeds by updating the stabiliser description of the state 
through the course of the computation and we get a description of the final 
state in addition to the output probabilities. We adopt here a different 
approach [3]. Suppose (wlog) that the final measurement is on the first line, 
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having outputs 0,1 with probabihties PQ,pi respectively. Then the difference 
Po — Pi is given by the expectation value of Zi = Z (g) 7 (g) . . . (g) 7 in the final 
state C \ipo): 

Po-Pi = (^o|CtZiC|^o) (2) 

This computation suffices to simulate the output (as also po + pi = 1). 
Now Zi is clearly a product of PauH operations so C^ZiC also has the 
product form Pi (g) . . . (g) P„ for Pauli operations Pi (whose identity can be 
determined in linear time by an update rule for successive conjugations by 
the elementary gates in the circuit). Hence if iV'o) = l^i) • • • l^n) is any 
product state we get 

n 
k=l 

which can clearly be calculated in time 0(n) (as a product of n terms of 
fixed size) giving an efficient (linear time) simulation of the Clifford circuit. 

The essential ingredients of the above proof are the following. 
(SIMl): we have a set Sn of n-qubit operations such that (V'ol 5'|V'o) can 
be computed in poly(n) time for any S E Sn and any allowed input state 

l^o); 

(For the GK theorem 5„ is the n-qubit Pauli group Vn-) 

(SIM2): we have a class /C„ of unitary operations such that K^SK G Sn 

for all S G Sn and K £ K-n- 

(For the GK theorem /C„ is the Clifford group C„.) 

Then if Zi is in 5„ for all n (or can be expressed in suitably simple terms 
using elements of 5„, c.f. later) then it follows (just as in the above outlined 
proof) that circuits of gates from /C„, with input state iV'o) and output 
measurement of Z on the first line, can be classically efficiently simulated. 

Note that this simulation result, resting on (SIMl) and (SIM2) does not 
actually require any special group (or other algebraic) structure on Sn or 
/C„. For example, the fact that Vn is a group is not needed at all in our proof 
of the GK theorem in contrast to the usual proof resting on the stabiliser 
formalism, depending heavily on the subgroup structure of Vn- 

Turning now to matchgates we will show that Valiant's theorem can be 
understood as just another example of the above formalism with a suitably 
clever choice of Sn and /C„. For n qubits we introduce the 2n Pauli product 
operators (omitting tensor product symbols (g) throughout): 

ci=X7...7 c^ = ZXI...I ■■■ C2k~i = Z ...Z XI ...I 

C2 = YI...I C4 = ZYI...I ■■■ C2k =Z...ZYI...I ^' 
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where X and Y are in the k slot for C2k-i and C2fc, and k ranges from 1 to n. 
For Sn we take the linear span of ci , . . . , C2n which is a 2n-dimensional vector 
space (in contrast to the group Vn)- Since each Cj is a product operator and 
a general vector v E 5„ is a linear combination of only 2n of them, it is 
clear that (V'ol iV'o) is poly(n)-time computable if IV'o) is a product state 
i.e. (SIMl) is satisfied. 

Next we can verify by straightforward direct calculation that if U is any 
n.n. G{A, B) gate then WcjU E Sn for all j so U^vU E 5„ for any ti E 5„ 
i.e. property (SIM2) is satisfied. More explicitly note that if [/ is a n.n. 
G{A,B) gate, it applies to two consecutive qubit lines so (from eq. ^) the 
part of Cj that it "sees" can only be one of 

ai = ZZ a2 = ZX 03 = ZY 04 = XI = YI or = //. (5) 

Then a straightforward calculation with 4 by 4 matrices shows that for each 
i, G{A, B)^ aiG{A, B) always returns a linear combination of allowable a^'s 
and property (SIM2) follows immediately. 

It is instructive to note that if we attempt to apply a G{A, B) gate on 
not nearest-neighbour qubit lines then in addition to the six terms in eq. ([5]) 
we can get a further possibility, namely = ZI on the chosen two lines. 
But now we can check that G{A, B)"^ ajG{A, B) does not generally lie in the 
span of the allowed Pauli products at those lines, and property (SIM2) is 
violated. This give a way of understanding the curious n.n. requirement for 
G{A^ B) actions in theorem [21 which has no analogue in the GK theorem 
(as Vn is defined by a uniformly local product requirement). 

With properties (SIMl) and (SIM2) we can say that if M is the total 
operation of any n.n. matchgate circuit on n lines then (V'ol Af^DM |^o) 
is poly(n)-time computable for any D ^ Sn- To complete our simulation 
theorem we want to set D = Z^ = I . . . I ZI . . . I (i.e. Z on the k^^ line) 
to obtain po — pi for a measurement on the k^^ line. In the GK theorem 
with Sn = Vn we had Z^ E Vn directly. In the present case we do not 
have Zk E 5„ but looking at eq. ([!]) we see that Zi = —iciC2 and generally 
Zk = —ic2k-iC2k- Then, for example, 

M^ZiM = -iM^ciC2M = -i{M^ ciM){M^ C2M) (6) 

and each bracket in the last expression is a linear combination of Cj's. Thus 
Po - Pi = (V'ol Af^-Z'iM IV'o) has the form X^jj (V'ol QCj iV'o)- Since 
the Ci's are product operators, so are the 0{v?) product terms CjCj in the 
final sum. Hence pq —pi is again poly(n)-time computable but now we have 
0{v?) terms instead of the previous 0{n) terms in the sum. This completes 
a proof outline of Valiant's theorem [2l 
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3 Extensions of simulatable circuits 



We may now view Clifford circuits and matchgate circuits as two "islands" 
of quantum processes in the space of all quantum computations, that offer 
no computational time benefit over classical computations. As such, it is 
interesting to try to characterise their relationship to the whole and one 
approach is to consider what (minimal) extra ingredient suffices to expand 
their computational power to regain full universal efficient quantum compu- 
tation. 

In the case of Clifford circuits it is well known (e.g. see pQ) that the 
inclusion of the phase gate \fP = diag(l, e^'^/^) suffices, and more generally, 
(using a result of Shi [10], noting that CNOT is a Clifford operation), the 
inclusion of essentially any single extra non-trivial 1-qubit gate will suffice. 

For the case of matchgate circuits we have the following intriguing result. 

Theorem 3 Let Cn he any uniform family of quantum circuits with output 
given by a Z basis measurement on the first line. Then Cn may be simulated 
by a circuit of G{A, B) gates acting on n.n. or next n.n. lines only (i.e. on 
line pairs at most distance 2 apart) with at most a constant factor increase 
in the size of the circuit. 

A proof of this theorem may be found in [2] and here we just make a few 
remarks. Comparing theorems [2] and [3] we see that the gap between classical 
and full quantum computational power can be bridged by a very modest 
use of a seemingly innocuous resource viz. the ability of matchgates to 
act on next n.n. - instead of just n.n. - qubit lines. Equivalently this 
may be characterised by use of the SWAP operation (on n.n. lines) in a 
very constrained context where ladders of consecutive SWAPs (which would 
allow 2-qubit gates to act on arbitrarily distant lines) are not even allowed. 
From this perspective, the power of quantum (over classical) computation 
is attributable to the mere inclusion of such isolated single SWAP gates. 
The result becomes perhaps even more striking if we note that SWAP 
itself is very close to being expressible in the allowed G{A, B) form. Indeed 
SWAP = G{I, X) and fails only through a mere minus sign in det X = 
— detJ. Thus if we drop the det^ = det-B condition in eq. ([T]), then the 
resulting G{A, B) gates acting on n.n. lines become efficiently universal for 
quantum computation. 

Is it conceivable that the passage from n.n. to next-n.n. use of G{A, B) 
gates may be achieved while maintaining classical simulatability? We may 
argue on formal complexity theoretic grounds that this is highly implausible. 
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Indeed it is shown in [2] that the classical complexity classes NP and PP 
(cf. [7]) would then become classically poly-time decideable i.e. we would 
get P=NP=PP (as well as P=BQP). Thus an extra supra-classical compu- 
tational power must be associated to the single distance extension of the 
range of n.n. 2-qubit G{A, B) gates in general matchgate circuits, if these 
classical computational complexity classes are to be unequal. 

4 Concluding remarks 

^From the viewpoint of (SIMl) and (SIM2) we see a formal similarity be- 
tween the GK theorem and Valiant's theorem although these results arose 
historically from very different considerations. This suggests that we might 
be able to construct further interesting classes of classically simulatable cir- 
cuits by simply taking other choices of Sn and identifying a suitable asso- 
ciated /C„. However "interesting" pairs {Sn-,K,n) appear to be difficult to 
invent - the known examples arising as outcomes of some prior elaborate 
underlying mathematical structures. In the GK case we have the identifi- 
cation of the Clifford group via a lengthy argument with group theoretic 
ingredients (see e.g. appendix in [11]) applied to the Pauli group Vn which 
is a well known structure in the subject. 

However in the case of Valiant's theorem, how might we initially come 
upon this result, and guess the choice for Sn that we used (i.e. eq. (|4]) 
and its linear span)!? Actually the operators in eq. ([!]) are well known in 
physics - they comprise the so-called Jordan- Wigner representation [i8j that 
appears in the theory of non-interacting fermions. The connection between 
Valiant's theorem and simulation of free fermions was recognised by Knill 
[5] and Terhal and DiVincenzo [6] and our proof of Valiant's theorem above 
is a re-writing of this connection. A more formal mathematical treatment 
(albeit without reference to fermions) based on abstract properties of the 
mathematical structure of Clifford algebras is given in [2] which also clar- 
ifies the appearance of matchgates as normalisers of the linear part of the 
Clifford algebra, leading to property (SIM2). We will not elaborate here on 
these further ingredients (detailed in ^) except to point out that again here, 
we have a significant underlying theory leading to the choice of 5„ and the 
identification of its associated normalisers /C„. Perhaps an intuitive signal 
feature of such an underlying theory is some construction that could po- 
tentially produce an exponentially large structure but surprisingly remains 
only polynomially complex. In the case of the Pauli group conjugation 
by arbitrary V £ C/(2") can generate general n-qubit matrices for which 
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the calculation of the expectation value in eq. ([2]) becomes exponentially 
inefficient. But the special case of V being Clifford guarantees a polynomial 
simplicity via the preserved product structure. In the case of the Cj's of eq. 
([U, conjugation by an arbitrary V S [7(2") leads to a general element of the 
full Clifford algebra generated by the the Cj's [2] - a space of exponential 
dimension 2^" - but again the special case of n.n. matchgates (associated 
to a theory of quadratic hamiltonians [2]) guarantees that the conjugates 
remain in the polynomially small subspace of linear elements of the full Clif- 
ford algebra. It is an interesting open problem to exhibit further examples of 
such simplifications and of our formalism (SIMl), (SIM2), that may already 
exist within the literature of the theory of some yet more general kind of 
algebraic structure. 
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